Detection of a flow induced magnetic field eigenmode in the Riga dynamo facility 
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In an experiment at the Riga sodium dynamo facility, a 
slowly growing magnetic field eigenmode has been detected 
over a period of about 15 seconds. For a slightly decreased 
propeller rotation rate, additional measurements showed a 
slow decay of this mode. The measured results correspond 
satisfactory with numerical predictions for the growth rates 
and frequencies. 

PACS numbers: 47.65. +a, 52.65.Kj, 91.25.Cw 



Magnetic fields of cosmic bodies, such as the Earth, 
most of the planets, stars and even galaxies are believed 
to be generated by the dynamo effect in moving elec- 
trically conducting fluids. Whereas technical dynamos 
consist of a number of well-separated electrically con- 
ducting parts, a cosmic dynamo operates, without any 
ferromagnetism, in a nearly homogeneous medium (for 
an overview see, e.g., (|] and 

The governing equation for the magnetic field B in an 
electrically conducting fluid with conductivity a and the 
velocity v is the so-called induction equation 



(9B 1 

— = curl(v x B) H AB 

at ixqij 



(1) 



which follows from Maxwell equations and Ohms law. 
The obvious solution B = of this equation may be- 
come unstable for some critical value Rm c of the mag- 
netic Reynolds number 



Rm = fioaLv 



(2) 



if the velocity field fulfills some additional conditions. 
Here L is a typical length scale, and v a typical velocity 
scale of the fluid system. Rm c depends strongly on the 
flow topology and the helicity of the velocity field. For 
self-excitation of a magnetic field it has to be at least 
greater than one. For typical dynamos as the Earth outer 
core, Rm is supposed to be of the order of 100. 

The last decades have seen an enormous progress of 
dynamo theory which deals, in its kinematic version, 



with the induction equation exclusively or, in its full ver- 
sion, with the coupled system of induction equation and 
Navier-Stokes equation for the fluid motion. Numerically, 
this coupled system of equations has been treated for a 
number of more or less realistic models of cosmic bodies 
(for an impressive simulation, see 

Quite contrary to the success of dynamo theory, ex- 
perimental dynamo science is still in its infancy. This 
is mainly due to the large dimensions of the length scale 
and/or the velocity scale which are necessary for dynamo 
action to occur. Considering the conductivity of sodium 
as one of the best liquid conductors (a w 10 7 (17m) -1 
at 100°C) one gets hqcf ~ 10 s/m 2 . For a very efficient 
dynamo with a supposed Rm c — 10 this would amount 
to a necessary product Lv = 1 m 2 /s which is very large 
for a laboratory facility, even more if one takes into ac- 
count the technical problems with handling sodium. His- 
torically notable for experimental dynamo science is the 
experiment of Lowes and Wilkinson where two ferromag- 
netic metallic rods were rotated in a block at rest [0. 
A first liquid metal dynamo experiment quite similar to 
the present one was undertaken by some of the authors in 
1987. Although this experiment had to be stopped (for 
reasons of mechanical stability) before dynamo action oc- 
curred the extrapolation of the amplification factor of an 
applied magnetic field gave indication for the possibility 
of magnetic field self-excitation at higher pump rates f| . 
Today, there are several groups working on liquid metal 
dynamo experiments. For a summary we refer to the 
workshop " Laboratory Experiments on Dynamo Action" 
held in Riga in summer 1998 S. 

After years of preparation and careful velocity profile 
optimization on water models, first experiments at the 
Riga sodium facility were carried out during November 
6-11, 1999. The present paper comprises only the most 
important results of these experiments, one of them being 
the observation of a dynamo eigenmode slowly growing 
in time at the maximum rotation rate of the propeller. 
A more comprehensive analysis of all measured data will 
be published elsewhere. 
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D 1=0.25 m 
D2=0.43 m 
D3=0.80 m 
L=2.91 m 



FIG. 1. The Riga dynamo facility. Main parts comprising: 
1 - Two motors (55 kW each), 2 - Propeller, 3 - Helical flow 
region, 4 - Back-flow region, 5 - Sodium at rest, 6 - Sodium 
storage tanks, * - Position of the flux-gate sensor, x - Posi- 
tions of the six Hall sensors. 

The principal design of the dynamo facility, together 
with some of the most important dimensions, is shown 
in Fig. 1. The main part of the facility consists of a 
spiral flow of liquid sodium in an innermost tube (with a 
velocity up to the order of 15 m/s) with a coaxial back- 
flow region and a region with sodium at rest surrounding 
it. The total amount of sodium is 2 m 3 . The sodium 
flow up to 0.6 m 3 /s is produced by a specially designed 
propeller which is driven by two 55 kW motors. 

All three sodium volumes play an important role in 
the magnetic field generation process. The spiral flow 
within the immobile sodium region amplifies the mag- 



The back-flow is 
The result is an 



nctic field by stretching field lines 
responsible for a positive feedback 
axially non-symmetric field (in a symmetric flow geome- 
try!) slowly rotating around the vertical axis. Hence, a 
low frequency AC magnetic field is expected for this con- 
figuration. Concerning the azimuthal dependence of the 
magnetic field which includes terms of the type exp(?m</?) 
with in general arbitrary m it is well-known that for 
those Rm available in this experiment only the mode with 
m = 1 can play any role S. A lot of details concerning 



the solution of the induction equation for the chosen ex- 
perimental geometry and the optimization of the whole 
facility in general and of the shape of the velocity profiles 
in particular can be found in Q, 0, and jicj ]. 

For the magnetic field measurements we used two dif- 
ferent types of sensors. Inside the dynamo, close to the 
innermost wall and at a height of 1/3 of the total length 
from above, a flux-gate sensor was positioned. Addition- 
ally, 8 Hall sensors were positioned outside the facility 
at a distance of 10 cm from the thermal isolation. Of 
those, 6 were arranged parallel to the dynamo axis with 
a relative distance of 50 cm, starting with 35 cm from the 
upper frame. Two sensors were additionally arranged at 
different angles. 

After heating up the sodium to 300° C and pumping it 
slowly through the facility for 24 hours (to ensure good 
electrical contact of sodium with the stainlcss-stccl walls) 
various experiments at 250° C and around 205° C at dif- 
ferent rotation rates of the propeller were carried out. 

According to our numerical predictions, self-excitation 
was hardly to be expected much above a temperature 
of 200° C since the electrical conductivity of sodium de- 
creases significantly with increasing temperature. Nev- 
ertheless, we started experiments at 250°C in order to 
get useful information for the later dynamo behaviour 
at lower temperature, i.e. at higher Rm. Although 
the experiment was intended to show self-excitation of 
a magnetic field without any noticeable starting magnetic 
field, kick-field coils fed by a 3-phase current of variable 
low frequency were wound around the module in order 
to measure the sub-critical amplification of the applied 
magnetic field by the dynamo. This measurement phi- 
losophy was quite similar to that of the 1987 experiment 
pd| and is based on generation theory for prolongated 
flows as the length of our spiral flow exceeds its diame- 
ter more then ten times. Generation in such a geometry 
should start as exponentially high amplification of some 
kick-field (known as convective generation) and should 
transform, at some higher flowrate, into self-excitation 
without any external kick-field. 

As a typical example of a lot of measured field ampli- 
fication curves, Fig. 2 shows the inverse relation of the 
measured magnetic field to the current in the excitation 
coils for a frequency of 1 Hz and a temperature of 205°C 
versus the rotation rate of the propeller. The two curves 
(squares and crosses) correspond to two different settings 
of the 3-phase current in the kick-field coils with respect 
to the propeller rotation. An increasing amplification of 
the kick-field can be clearly observed in both curves until 
a rotation rate of about 1500 rpm. These parts of both 
curves point to about 1700 rpm which might be inter- 
preted as the onset of convective generation [|| , j|] . If the 
excitation frequency would be exactly the one the system 
likes to generate as its eigenmode, the curves should fur- 
ther approach the abscissa axis up to the self-excitation 
point. As 1 Hz does not meet exactly the eigenmode fre- 
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quency the points are repelled from the axis for further 
increasing rotation rates as it is usual for externally ex- 
cited linear systems passing the point of resonance (for 
this interpretation, see also Fig. 5). 
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FIG. 2. Dependence of the magnetic field amplification on 
the propeller rotation rate for T=205°C and f=l Hz. The 
ordinate axis shows the inverse relation of the measured mag- 
netic field to the current in the kick-field coils. Squares and 
crosses correspond to two different settings of the 3-phase 
current in the kick-field coils with respect to the propeller 
rotation. 

It should be underlined that all points on Fig. 2 ex- 
cept the rightmost one are calculated from sinusoidal field 
records showing the same 1 Hz frequency as the kick-held. 
However, the rightmost point at 2150 rpm is exceptional. 
Let us analyse the magnetic field signal at this rotation 
rate in more detail. Fig 3a shows the magnetic field mea- 
sured every 10 ms at the inner sensor in an interval of 15 
s. Evidently, there is a superposition of two signals. 

Numerically, this signal (comprising 1500 data points) 
has been analyzed by means of a non-linear least square 
fit with 8 free parameters according to 

B(t) = A ie pit sin (27r/it + fa) + A 2 e P2t sin (27r/ 2 f + fa) 

The curve according to this ansatz (which is also shown 
in Fig. 3a) fits extremely well into the data giving the 
following parameters (the errors are with respect to a 
68.3 per cent confidence interval): 

Ai= (0.476 ± 0.004) mT, Pl = (-0.0012 ± 0.0003) s" 1 
h = (0.995 ± 0.00005) s -1 , fa = -0.879 ± 0.012 

A 2 = (0.133 ± 0.001) mT, p 2 = (0.0315 ± 0.0009) s" 1 
h = (1-326 ± 0.00015) s" 1 , fa = (0.479 ± 0.009) 

The positive parameter p 2 = 0.0315 s _1 together with 
the very small error gives clear evidence for the appear- 
ance of a self-exciting mode at the rotation rate of 2150 
rpm. Fig. 3b shows in a decomposed form the two con- 
tributing modes, the larger one reflecting the amplified 



field of the coils and the smaller one reflecting the self- 
excited mode. 
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FIG. 3. Measured magnetic field and fitting curve (a). De- 
composition of the fitting curve into two curves with different 
frequencies (b). 

For reasons of some technical problems, this highest ro- 
tation rate could be hold only for 20 seconds after which 
it fell down to 1980 rpm. At that lower rotation rate the 
coil current was switched of suddenly. Figure 4 shows 
the magnetic field behaviour at three selected Hall sen- 
sors positioned outside the dynamo. This mode has a fre- 
quency of / = 1.1 s _1 and a decay rate of p = —0.3 s _1 . 
A similar signal was recorded by the inner fluxgate sen- 
sor, too. 

It is interesting to compare the frequencies and growth 
or decay rates at the two different rotation rates 2150 rpm 
and 1980 rpm with the numerical predictions. These are 
based on the outcomes of a two-dimensional time depen- 
dent code which was described in As input velocity 
for the computations an extrapolated velocity field based 
on measurements in water at two different heights and at 
three different rotation rates (1000, 1600, and 2000 rpm) 
was used. Fig. 5 shows the predicted growth rates and 
frequencies for the three temperatures 150°C, 200°C, and 
250°C which are different due to the dependence of the 
electrical conductivity on temperature. The two pairs of 
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points in Fig. 5 represent the respective measured val- 
ues. Having in mind the limitations and approximations 
of the numerical prognosis || the agreement between pre- 
calculations and measured values is good, particularly re- 
garding the frequencies of the magnetic field cigenmodc. 
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FIG. 4. Magnetic fields measured at 3 selected positions 
outside the dynamo module after switching off the coil cur- 
rent. 



the observed growth rate was still very small. The corre- 
spondence of the measured growth rates and frequencies 
with the numerical prognoses is convincing. The general 
concept of the experiment together with the fine-tuning 
of the velocity profiles || have been proven as feasible 
and correct. The facility has the potential to exceed the 
threshold of magnetic field self-excitation by some 20 per 
cent with respect to the critical magnetic Reynolds num- 
ber. The experiment will be repeated at lower temper- 
ature when the technical problems with the seal will be 
resolved. For lower temperature, a higher growth rate 
will drive the magnetic field to higher values where the 
back-reaction of the Lorentz forces on the velocity should 
lead to saturation effects. 
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The main part of the experiment was originally 
planned at T=150°C where self-excitation with a much 
higher growth rate was expected. Unfortunately, the 
safety rules required to stop the experiment at T=205°C 
since technical problems with the seal of the propeller 
axis against the sodium flow-out have been detected. It 
is worth to be noted that the overall system worked stable 
and without problems over a period of about five days. 
The sealing problem needs inspection, but represents no 
principle problem. 
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FIG. 5. Numerical predictions for growth rates p and fre- 
quencies / of the dynamo eigenmode in dependence on the 
rotation rate for three different temperatures, and measured 
values. 

For the first time, magnetic field self-excitation was ob- 
served in a liquid metal dynamo experiment. Expectedly, 
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